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Abstract 

We discuss the generation of parametric X-rays in the new photoinjector at the 
FAST (Fermilab Accelerator Science and Technology) facility in Fermilab. These ex¬ 
periments will be conducted in addition to channeling X-ray radiation experiments. 
The low emittance electron beam makes this facility a promising source for creating 
brilliant X-rays. We discuss the theoretical model and present detailed calculations of 
the intensity spectrum, energy and angular widths and spectral brilliance under differ¬ 
ent conditions. We also report on expected results with parametric X-rays generated 
while under channeling conditions. 


1 Introduction 

Energetic charged particles traveling through a crystal can produce X-rays by several mech¬ 
anisms. Incoherent bremsstrahlung and transition radiation give rise to a continuous spec¬ 
trum while channeling radiation (CR) and parametric X-ray (PXR) radiation produce quasi- 
monochromatic discrete X-ray spectra. One of the main advantages of using crystals is that 
CR and PXR produce hard X-rays with much lower energy electrons compared to, for ex¬ 
ample, synchrotron radiation produced X-rays in circular rings, . It takes a 3 GeV electron 
beam (assuming a bend field of 1 T) to generate X-rays with a critical energy of lOkeV 
via synchrotron radiation while 10 MeV electrons have sufficed with channeling and para¬ 
metric radiation at the same energy. Hard X-ray generation using crystals and 50 MeV 
electrons is one of the planned set of experiments at Fermilab’s L-band photoinjector in the 
FAST facility (formerly called ASTA) [[IlEl, currently being commissioned . The major 
goal of these experiments is to demonstrate that such a photoinjector with a low emittance 
electron beam can serve as a model for a brilliant compact X-ray source when scaled to a 
higher gradient X-band photoinjector. 

The detailed characteristics of CR expected at FAST was discussed in Ol. In this pa¬ 
per we will consider the spectral brilliance of PXR under various conditions at FAST. The 
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PXR mechanism was first discussed several decades ago (411511, experimentally measured 
first with electrons in 1985 O and since then observed at many laboratories; several re¬ 
views are now available (TUlll. PXR has also been observed from 400 GeV protons using 
a bent crystal at CERN’s SPS accelerator (91. The characteristics of PXR differ from CR 
in several ways. In CR emission, the X-ray energy spectra is discrete at electron ener¬ 
gies below 100 MeV and the frequencies depend on the particle energy, while in PXR 
they are independent of particle energy at relativistic particle speeds. PXR can be gener¬ 
ated at large angles from the particle’s direction which differentiates it from both CR and 
bremsstrahlung making the background contribution significantly less than the signal. The 
PXR spectral lines are also more monochromatic than CR, the width is at least an order 
of magnitude smaller. The disadvantage of PXR is that the photon yield is about two to 
three orders of magnitude smaller than that of CR. On the other hand, PXR can be gen¬ 
erated simultaneously with CR generation thus potentially allowing multiple X-ray beams 
with different spectra and in different directions. The crystal requirements for PXR and CR 
production are similar, namely high thermal conductivity, low photon absorption length, 
high dielectric susceptibility and large lattice spacing. In this article we will consider the 
spectral brilliance from a PXR source under different conditions. In general, the brilliance 
will be several orders of magnitude lower than that from the brightest X-ray sources such as 
XFELs or inverse Compton scattering. However, compared to those sources, a PXR source 
can deliver X-rays suitable for industrial and medical applications with significantly lower 
cost, complexity and size. The special feature of the photoinjector at FAST for generating 
brilliant X-rays using crystals is that low emittances can be generated by shaping the laser 
spot size on the cathode and even lower emittances have been obtained with field emission 
(FE) nanotip cathodes (TTl . In addition, the photocathode and FE cathodes operate at GHz 
frequencies, so the high repetition rate allows for low bunch charge required for these low 
emittances. 

In Section 2, we discuss the PXR spectrum, notably the photon energy and the spectral 
distribution dependence on the crystal geometry. In Section 3, we present calculations of 
the energy width with contributions from geometrical effects and multiple Coulomb scat¬ 
tering, while in Section 4 we consider angular broadening and in Section 5 the spectral 
brilliance. Section 6 contains specific calculations for the photoinjector at FAST; this in¬ 
cludes cases with the use of the present goniometer holding the crystal, the options with 
a new goniometer, and finally PXR emission under channeling conditions. We conclude 
with a summary of results in Section 7. 


2 Characteristics of the PXR spectrum 

PXR emission occurs when the virtual photons accompanying the charged particle scatter 
off the atomic electrons in the crystal and interfere constructively along certain directions. 
The incoming virtual photon’s wave-vector k, and the outgoing real photon’s wave vector 
ky are related by the Bragg condition for momentum transfer 

ky = k,-hmg (1) 
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where g is the reeiproeal lattiee veetor of the seattering planes and m is an integer. Writing 
k, = (0j8/v, \if = (oCl/ (c/ y/e), and taking the sealar produet of the above equation with jS 
yields the outgoing PXR photon energy as 


E 


h(D 




( 2 ) 


Here v = j8/3c is the veloeity vector of the particle, is the unit velocity vector. For 

A 

electron energies in the range of tens of MeV, we can approximate p ~ 1. is the unit 
vector along the direction of the emitted photon and e is the real part of the permittivity. If 
the first order spectrum is obtained by reflection from a plane with Miller indices (h,k,/), 
multiples of this frequency occur from reflections off planes with indices {mh, mk, ml) with 
m > 1. The above equation can also be derived by requiring that the phase difference 
between photons reflected from adjacent lattice planes be an integer multiple of 2;r. At 
X-ray energies, the frequency dependent real part of the dielectric function can be written 
as £{(o) ~ 1 — {(0p/(0)^ The plasma frequency (Op for most crystals is in the range 10-90 
eV while X rays have keV range energies. Approximating e 1, the first order (m = 1) 
PXR energy can be written as 
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where 6b is the angle of the crystal plane with the beam direction, a is the angle of the 
electron with the central electron beam direction, and 6o is the observation angle of the 
emitted radiation with the beam direction. This shows that the photon energy is indepen¬ 
dent of beam energy and can be changed by rotating the crystal with respect to the beam 
direction. 

The expression for the angular intensity distribution calculated from a kinematic theory 
can be written as lH^ [T3l [T4ll 
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Here Ox, By are the angular deviations from the direction of specular reflection (6d = 2Ob) 
in the planes parallel and perpendicular respectively to the diffraction plane, OCy is the 
fine structure constant. La is the photon absorption length in the crystal, t is the crystal 
thickness, is the frequency dependent part of the Fourier transform of the dielectric 
susceptibility and depends on the plasma frequency (Op, the atomic number Z of the crystal 
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and on 5(g),F(g). 5(g) is the crystal structure function, the sum is evaluated over the atom 
location r in a unit cell. For FCC cubic crystals like diamond and silicon, 5(g) 7 ^ 0 only if 
the Miller indices (h,k, Z) of the reflection plane are either all odd or all even. If they are 
all even, must be divisible by four. For an amorphous material 5(g) = 0, hence 

there is no PXR emission. F’(g) is the atomic scattering form factor which is the Fourier 
transform of p, the density distribution of the atomic electrons. The Debye-Waller factor 
accounts for thermal vibrations and is close to unity at room temperature. 

Among the assumptions made in deriving this result are that the beam is relativistic 
(7 3> 1), deviations of the photons from the Bragg specular reflection angle are small 
( 0 x, 0 j< 1 ) and the crystal thickness t is larger than the photon formation length. 

The angular factor, the term in square brackets on the right hand side of Eq. ([4]) contains 
the dependence on the angles By, the beam energy 7 and on the X-ray frequency CO. The 
two terms in the numerator of the angular factor describe the contributions from the two 
polarizations. If the observation angle 20b = ^/2, the horizontal contribution vanishes 
and the observed spectrum is completely vertically polarized at all angles Ox^Oy When 
projected onto one of the two angles the angular distribution can be either single peaked 
or double peaked, depending on the value of the orthogonal angle. The extrema of the 

angular distribution in 9x are at 9x — 9^^^ + 9^{1 —Iscc^IBb) while those in 9y are 

at 0^ = Oj 9pf^ + 9x{l — 2cos^20g). When there is a single peak, the maximum is at 

zero, while with double peaks, the maxima are close to ±9ph and there is a local minimum 
at zero angle. 

The angular distribution depends on the beam energy only through the angle Bp^ which 
occurs in the angular factor. At low beam energies such that y (co/cOp), we have Bph ~ 
1/7 and the intensity increases as 7 ^. This can also be seen by a power series expansion of 
the angular factor in terms of the parameter yl^cojoop) 
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The denominator (1 + 7^(0^ + By) is of order unity over the useful range of angles 0^, By. 
As the beam energy increases, the angular factor and the angular intensity distribution 
reaches a maximum around 7 = cajcOph, levels off and then decreases slowly at higher 
beam energy. This behavior can be seen in Fig. [T|in which the value of the angular factor 
at 0;i; = 0, 0y = 1 /7 is plotted. In the lower of the two curves where 7 ranges from zero 
to larger than cojcap the intensity levels off, while in the upper curve, 7 < ( 0 /( 0 p over the 
entire range, so the intensity grows monotonically 

The atomic structure function E(g) and consequently the susceptibility are calculated 
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Figure 1: Dependence of the angular factor with the electron beam’s relativistic factor y 
for two values of co/cOp with cOp = 38eV. The lower value co/cOp corresponds to an X-ray 
energy of 5.7 keV and the higher value to a X-ray energy of 19keV. 


from the expressions 

Fis)=Ms) + if + if"), Ms) = ^a,exp[-h,-5^] -be, 5 = Ang 

i 

Xo = + if], Xg = -^^[(Ms) + / - Z) + if"]Shki (9) 

The coefficients {at^bi^c) are the Cromer-Mann coefficients lITSl while the frequency de¬ 
pendent form factors {f ,f ) can be obtained from a database maintained by NIST IfT^ . 
Here rg is the classical electron radius, Nc is the number of atoms in the unit cell, Vc is the 
volume of the unit cell, Z is the atomic number of the crystal and Shki is the crystal structure 
factor for the plane with indices (/?,k, /). The photon attenuation length at a wavelength A 
can be found from = X/{2K\lm{XQ)\)■ 

We now discuss the geometric factor fgeo in Eq.dH). The unit vectors fi, v, Cl are defined 
and shown in FiglU The figure on the left in Figure [2] shows the electron velocity vector 
v, the normal n to the crystal surface, and the normal g to the crystal planes. The figure 
on the right in Figure [2] shows the direction of the detector’s central axis Cl, and the angles 
between the vectors. It is clear that in order for the crystal plane to be a reflecting plane, the 
angle nj2 — Ob between g and -v and also the angle ^ between g and fl must satisfy 0 < 
{n 12 -Ob AX) ^kI2 where C, is the angle between gand n. Bragg geometry corresponds 
to ^ = 0 while in Laue geometry ^ = ;r/2. If is the angle of the detector relative to the 
velocity vector v, then we have 

(j) = ^-{0 d-9b), n-v =-sin(0B-t-C), n-Cl = sm{0D - 9 b - Q 

The geometric factor fgeo in the intensity expression can then be written as 
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Figure 2: Left: A general erystal geometry for an arbitrary orientation of the refleeting 
plane. Reeiprocal lattiee veetor g: normal to the erystal plane, n: normal to the erystal 
surfaee. Bragg geometry corresponds to the crystal plane which has g||n, Laue geometry 
to g ± n. Right: The different angles between g, n, v: electron velocity vector, and Q: 
direction of the detector’s central axis 

Note that if 8 d = 8b + the geometric factor fgeo = 0 because the photons travel along 
the larger transverse dimensions of the crystal and will be mostly absorbed in the crystal. 

We mention here that a refinement to the kinematical theory is the dynamical diffraction 
theory ifTTl [TSl [T9ll which takes into account the coupling between the photon fields with 
wave vectors ki and kf + g via interaction with the crystal. This coupling gives rise to 
additional PXR photons emitted in the forward direction in close proximity to the electron 
beam. This forward PXR was observed in experiments at the Mainz laboratory [l20l but care 
was required to extract this PXR emission signal as transition radiation and bremsstrahlung 
are also emitted in the same direction. We will not discuss this forward PXR emission here 
as it does not offer the relatively background free property of PXR emission at the Bragg 
angle. 


3 Energy spectrum broadening 

The intensity spectrum given by Eq.© predicts a delta function spectrum at integer multi¬ 
ples mcog of the Bragg frequency. In practice, there are several mechanisms which broaden 
the frequency of each line in the spectral distribution. We discuss the important sources 
and present analytical results for their contributions to the energy width and compare them 
to previous experimental results. 

3.1 Geometric effects 

From Eq.® , it follows that the photon energy depends on the angle of the incident electron 
and the photon direction. The incident angle will have a finite spread due to the beam 
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Figure 3: Effect of the beam spot size (a) and crystal thickness (b) on the angular photon 
spread. In (a), the beam spot diameter is represented by sq and sd is its projection on the 
detector. In (b), the crystal thickness is t, the effective thickness traversed by the electrons 
is tgff and to the projection on the detector. 


divergence while the photon angle hitting the detector can have a spread due to several 
effects including the finite beam spot size on the crystal surface, the detector size and the 
crystal thickness. 


Writing the spread in electron incident angle in the diffraction plane as Actx and the 
photon angle as Qd = 29b + 9x, we have from Eq.© to first order in Aa^ 

1 




2 sin^ 0B 


■ drACCr 


( 11 ) 


We set AUx = <3^ the beam divergence in the diffraction plane and 9x — Ax^etlR where Ax^et 
is the width of the detector and R is the distance from the crystal to the detector. 

Again from Eq.©, it follows that the energy spread is related to the photon angular 
spread as 

A0d 


^ E ^ tanPg 
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The impact of the finite beam size on A0d can be seen in Eig©. A finite size on the crystal 
surface projects to a spot size on the reflecting planes from which photons with a spread 
of angles can reach the detector. Eet sq be the beam spot diameter on the crystal and Sp its 
projection on the crystal’s reflecting plane. The angle the velocity vector v makes with the 
crystal surface is hence it follows from the figure that 
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On reflection from the crystal plane, the size Sp projects to a size sq = SpSm{0D — 0 b)- The 
angular spread in the photon angles resulting from the beam size is 
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where we have replaeed sq by the rms beam size Ox. 

The impaet of the finite erystal thiekness is seen in Figj^?. The effeetive erystal thiek- 
ness teff projeets to a length on the deteetor to = teffsmQo = tsmdo/ sin(0B + Q where 
t is the crystal thickness. Hence the angular spread in photon angles due to the crystal 
thickness is 

\ t sinOz) 

(A0z)Wf«/- ^ “i?sin(0B + C) 

Finally the detector size results in an angular spread {^dD)det = ^det/R- 


Adding the independent sources of beam angle spread and photon angle spread in 
quadrature, the energy spread due to these geometric effects is 


^ Ard^)2 ^^ot2 I (a^sin(0s + + {t • + i^dei)^ \ 

2sm Ob { sm(yg + c,j J 

(15) 

Here we have set 9o = 20 b, the direction for specular reflection. For typical beam and 
crystal parameters, the dominant contributions are from the beam spot size and the de¬ 
tector width while the contributions from the beam divergence and crystal thickness are 
significantly smaller. 


AE _ 1 
~E~R 


3.2 Multiple Scattering 


The contribution of multiple Coulomb scattering can be analytically calculated from the 
differential angular spectrum. Using Eq.®, it follows that the differential angular spectrum 
per unit length is given by 
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The spread in angles 0x,0y changes as the electron beam propagates through the crystal due 
to multiple scattering. Assuming the multiple scattering process to be Gaussian, the angles 
are sampled from distributions 


/(^x) = , exp[ 
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and a similar expression for the distribution in Oy. Here (7^, o'y are the rms beam divergences 
which increase with z as the beam propagates through the crystal. Writing the initial beam 
divergences as (ct^ q, q), the z dependent divergences are 


= [{c>'x,Q? + {aMs{z)) 
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and similarly for <7^,. Here Oxms Is the rms multiple scattering angle in the x direction. This 
can be found from the expression [|2T|| 
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where z is the path length traversed in the erystalj^g is the eleetron beam energy in MeV 
and Xrad is the radiation length. 


The multiple seattering weighted distribution funetion in Ox is then 
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Here ineludes all factors which do not depend on z, 0x, Qy Within the integrand, E{9x + 
(px) denotes that the photon energy is evaluated at the angle ( 9 x + 0x)- From this weighted 
distribution, the average and rms width of the energy spectrum can be found as 


_jEdN _ jE(dN/d9x)Msd9x 
^ ' JdN JidN/d9x)Msd9x 

(19) 

We will use Eq. (fT9l) to estimate the energy width due to multiple scattering instead of a 
Monte-Carlo simulation that is often used. 


There is another contribution to the linewidth from the photo-absorption in the crystal. 
Assuming a point source electron beam and no imperfections or multiple scattering, the 
photon wave train emitted by the beam has an intrinsic energy width given by [l22l 
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where Xo,i Is the imaginary part of the mean dielectric susceptibility. The intrinsic width 
has a minimum value for backward emission when 9b = %I2. This contribution is typically 
orders of magnitude smaller than the other contributions discussed above. 


The expressions for the energy width have been checked against measured values from 
a couple of earlier experiments, one at a low beam energy of 6.8MeV ll23ll and the other 
with beam energy of 56 MeV [|2^ . close to the FAST beam energy of 50 MeV. Since 
multiple scattering is more important at lower energy, comparison with the low energy 
result are a good check of the width from multiple scattering while the second case will be 
a good check of the geometrical width. Table [U shows the results of the comparison. The 
geometrical linewidth was calculated using Eq.(fT5l) and the multiple scattering contribution 
using Eqs.tfTSl) and (fT9l) . The two were then added in quadrature to yield the theoretical 
value shown in Table [IJ 


4 Angular spectrum broadening 

The measured angular intensity distribution represents a convolution of the intrinsic PXR 
intensity with the Gaussian response of the detector angular resolution, the beam diver- 
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Crystal 

Geometrical 

[eV] 

Multiple Scattering 
[eV] 

Intrinsic 

[eV] 

Theory(total) 

[eV] 

Experiment 

[eV] 

C(lll) 

35.9 

34.8 

2.1x10-4 

49.4 

51 

Si(220) 

524.9 

61.3 

2.8x10-4 

529 

540 ± 120 

Si(400) 

97.1 

9.7 

2.1x10-3 

98 

134 ± 56 


Table 1: Comparison of theoretical energy widths with experimental values after removing 
the effects of the detector energy resolution. First row: Experiments with 6.8 MeV elec¬ 
trons and a diamond crystal [|^ . Second and third rows: Experiments with 50-60 MeV 
electrons and silicon crystals |[2^[25]l . These values are taken from Table 36 in ll24l . The 
theoretical estimates include geometrical, multiple scattering and intrinsic contributions 
added in quadrature. 


gence and multiple scattering. Hence the measured intensity is of the form 
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where Ac is a constant to ensure photon number conservation after the convolution, d^N/dOxOy 
is the angular distribution without convolution, and the total angular resolutions cj, oj are 
given by 





( 22 ) 


where u = (x,y), Od,u is the detector resolution, R is the distance of the detector from the 
crystal, o'^^ is the beam divergence at the crystal, is the effective multiple scattering 
angle averaged over the path length tj cos Qb- 


In the absence of the broadening due to the convolution, the intrinsic angular width of he 
PXR spectral distribution is given by AOpxr — Op^. In crystals with thicknesses comparable 
to the attenuation length, the broadening due to multiple scattering is significant and the 
characteristic double peaked angular distribution is replaced by a broadened single peak 
distribution with the center filled in. 


Other sources of angular broadening are the higher order reflections from planes with 
spacings which are integer sub-multiples of the primary plane. These higher order reflec¬ 
tions produce higher photon energies with lower yield leads to a broadened distribution 
when recorded on a detector which sums over all photon energies. 


5 Spectral Brilliance 


The spectral brilliance of the photon beam is defined as the number of photons emitted per 
second per unit area of the photon beam per unit solid angle per unit relative bandwidth 
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Expressed in eonventional light souree units, the average speetral brillianee ean be written 
in terms of the averaged beam parameters and differential angular intensity speetrum per 
eleetron in a 0.1% bandwidth 
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where lav is the average eleetron beam eurrent, Ey is the energy of the X-ray line and Oy 
is the X-ray beam spot size. In the seeond line dN/dQ. is the angular yield in units of 
photons/(el-sr), and we set the photon spot size to the eleetron beam spot size in the erystal, 
i.e. Oy = Cg. Here AEy includes only the contributions to the spectral width from the crystal 
but not that from the detector resolution. 


Due to multiple scattering within the crystal, the electron beam divergence, beam size 
and emittance will grow as the electrons move through the crystal. Writing £m as the 
normalized electron emittance, the emittance growth as a function of the distance traversed 
z is 

^£n,x{z) = (26) 

where Px is the optical beta function at the crystal in the x axis. If the initial beam size 
at the crystal is Oxfi and the normalized emittance at the crystal entrance is then 

Px = y(o^Q/£N,o- The beam divergence grows as o'^{z) = \J (o-g(0))2 {o^^{z)Y and the 

average of the inverse beam size squared follows from 


K ^ y_i_L_\ = T M . _1_ 

PxEn{z) PxteffJo ^ [£Nfi + rl^x{(yMs(^))^] 


(27) 


This averaged expression will be used in Eq.(i25]) for the average brilliance. If the initial 
beam divergence and beam emittance are small compared to their increase through the 
crystal, and neglecting the logarithmic correction to the multiple scattering angle , we have 
(Oe(z), Og(z) y/z and the emittance grows as £n{z) z. 

Since the brilliance scales linearly with the electron current but inversely as the square 
of the electron spot size, it is advantageous to geneate as small a spot size as feasible even 
at the cost of reducing the beam current. Since the emittance grows with crystal thickness 
faster than the yield does, the maximum brilliance will also require the use of thin crystals, 
as will be shown later with numerical examples. 


6 PXR in FAST 

In the EAST beamline, a goniometer on loan from the HZDR facility, described in [l26ll . is 
presently available. It has two ports through which radiation can be extracted - one along 
the beam axis which will be used for channeling radiation and another at 90° to the beam 
axis which can be used to extract PXR. This determines that with the detector angle at 
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Parameter 

Value 

Units 

Beam energy 

50 

MeV 

Bunch charge 

20 

pC 

Length of a macropulse 

1 

ms 

Number of bunches/macropulse 

2000 

- 

Macropulse repetition rate 

5 

Hz 

Bunch frequency 

2 

MHz 

Interval between bunches 

0.33 

/is 

Bunch length 

3 

ps 

Crystal, thickness 

Diamond, 168 

/im 


Table 2: Bunch and macropulse parameters in FAST 


Plane 

X-ray energy 
[keV] 

^a,C 

[cm] 

La,air 

[cm] 

Attenuation 
in air 

Yield 

[photons/el-sr] 

AE 

[eV] 

(1,1,1) 

4.26 

0.0097 

12.72 

3.8x10-^ 

3.7x10"^ 

59 

(2,2,0) 

6.95 

0.043 

57.2 

0.17 

9.9x10^5 

93 

(4,0,0) 

9.83 

0.120 

144.9 

0.50 

8.8x10^5 

131 


Table 3: Photon yields and linewidths at a Bragg angle of 45°, observation angle of 90° 
in FAST from PXR off some planes. Bragg geometry in all cases, crystal thickness = 
0.168mm. The relative energy width in all cases is about 1%. The yield value includes the 
effect of attenuation over aim long path in air from the crystal to the detector. 


90 degrees, the Bragg angle must be 45 degrees in order to generate PXR with sufficient 
intensity. 

The goniometer already has a diamond crystal inside with its surface cut parallel to the 
(1,1,0) plane to generate channeling radiation from this plane. It has a thickness of 168/im 
and this will be the assumed value for most calculations reported here. In order to limit 
heating the crystal by the beam, we will assume low current operation with an average 
beam current of 200nA and a bunch charge of 20pC. At such charge values, low transverse 
emittance of the order of 100 nm can be obtained by suitably shaping the laser spot on the 
cathode If27l or alternatively with field emission cathodes [fTTlI . The main parameters of 
FAST and the crystal are shown in Tabled We choose R = 1 m as the crystal to detector 
distance and the active size of the detector plate to be 2cm x 2cm. 

Table [3] shows the PXR photon energies, the angular yield and the energy width from 
reflection off three of the possible low order planes with a 50 MeV electron beam. The 
yields include the effect of attenuation in air from the crystal to the detector. The recip¬ 
rocal lattice spacing g between adjacent planes with indices {h,k,l) is found from g = 
{Ik/ a)\/h^ + where a is the length of a unit cell. Consequently both the energy and 

absolute linewidths increase with increasing order. The yields are higher for the (2,2,0) and 
the (4,0,0) planes primarily due to the higher susceptibility Xg- Table |4] shows linewidth 
contributions from geometric effects and multiple scattering for each of the planes. The 
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Plane 

Geometrical Width (eV) 

Multiple-Scattering(eV) 

Total (eV) 

(111) 

42.8 

39.3 

59.1 

(220) 

69.5 

62.0 

93.2 

(400) 

98.3 

86.7 

131.1 


Table 4: Contributions to the energy width from geometric effects and multiple scattering. 


0 10 


0 05 


- 0.05 


-010 

-0 10 - 0.05 0.00 0.05 0 10 -0 10 - 0.05 0.00 0.05 0 10 

Figure 4: Contour plots of the differential yields as functions of (Ox, By) without broadening 
(left) and with broadening(right). The white areas in the left figure correspond to the highest 
intensities. The convolution lowers the peak intensities and results in a central broad peak 
instead of the two peaks without convolution. Both figures are for the (2, 2, 0) plane with 
a Bragg angle of 45° and a beam divergence of Imrad. These contour plots correspond to 
the images that would be observed on a detector plate. 

two effects are comparable in these cases. 

Figure |4] shows the two dimensional contour plots of the angular intensity spectrum 
projected on the (Ox, Oy) axes without and with convolution. The broadening effects are 
clearly visible as the two distinct peaks in the left figure merge into a single wider maxima 
in the right figure. 

For each primary plane, PXR emission also occurs from higher order planes at higher 
energies and lower intensities. For the (1,1,1) plane, the next allowed higher order plane is 
the (3,3,3) plane, since S(g) =0 for the (2,2,2) plane. The photon energy from the (3,3,3) 
plane is 12.8 keV with significantly reduced attenuation in air and resulting in a angular 
yield at the detector about two orders of magnitude higher than that from the (1,1,1) plane. 
For the (220) plane, second order reflections from (440) are allowed with photon energy of 
13.9keV. The (440) plane has an angular yield about a third smaller than the (220) plane 
even after including smaller attenuation at the higher energy. For the (400) plane, the 
second order reflection from the (800) plane produces 19.7keV photons and an angular 
yield about 10% that of the first order yield. The broadening of the angular distributions 
from the higher order reflections in these cases is found to be insignificant. 
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Crystal thickness/Attenuation length Crystal thickness/Attenuation length 


Figure 5: Left: Yield per eleetron per steradian as a funetion of the erystal thiekness relative 
to the photon attenuation length in the crystal for two different planes: (220) and (400). 
Photon attenuation lengths in the crystal are: 431 jum at 6.9keV, and 1198/im at 9.8keV. 
Right: Spectral brilliance as a function of the relative crystal thickness for the two planes. 


Figure [5] shows the impact of the crystal thickness on the angular yield and spectral 
brilliance for two different planes. The thickness is shown relative to the photon attenuation 
length, this length is larger for the (400) plane because of the higher photon energy. For the 
same absolute thickness, the angular yield with the (220) plane is larger, as seen in Table [3] 
for a single thickness but the angular yield for the same relative crystal thickness is higher 
with the (400) plane because the absolute thickness is larger. In both cases, the angular 
yield appears to saturate at a thickness of about 1.2La. The spectral brilliance for the same 
relative crystal thickness is larger with the (2,2,0) plane because the average emittance over 
the crystal is smaller with a smaller absolute thickness. In both cases, the brilliance reaches 
a maximum around (0.2 - 03)La. These plots show that the optimum crystal thickness 
depends on whether the photon yield or the spectral brilliance is the object of interest. 

6.1 New Goniometer 

The initial set of experiments will be conducted with the goniometer described above with 
the two ports. There is another goniometer under construction which will have a total of 
five ports through which X-rays could be extracted. A schematic of this new goniometer 
is shown in Figure 0 These ports will offer the opportunity of generating PXR at differ¬ 
ent Bragg angles and hence different energies, yields and spectral brilliance from those 
considered in the previous subsection. 

The angle of the port determines the direction of the photon emission and hence the vec¬ 
tor n. We write the unit vector along the detector axis as ^ = (tOi, (O 2 , -I- col + • 

If the PXR plane has Miller indices {h,k,l), the unit normal g can be one of ±{h,k,l)/ y/h^ + k^ + 1^. 
The conditions for reflection require that both g ■ (—v) > 0 and g ■ > 0. The first of these 

results in choosing the sign of g while the second determines if a chosen detector angle is 
suitable for observing reflected photons. Together the two conditions ensure that the in¬ 
cident beam and the reflected photons are on the same side of the reflecting plane. For a 
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7.9- 18.2 keV 



3.3 - 7.7 keV 


7.1 - 16.4 keV 


j 4.3 - 9.8 keV 


electrons' 3.2-7.5keV 


Figure 6: Schematic of the new goniometer with the five labeled ports. The X-ray energies 
expected from these ports with different planes, see Tabled are also shown. Note: Ports 
3 and 5 are not in the horizontal plane but are 45° to the horizontal. There is an exit port 
under port 5 and in the horizontal plane but not shown here. 

given plane with a normal g, the Bragg angle is determined from Ob = arcsin[—v-g]. For 
the goniometer ports shown in Figured the angle vectors to the detector locations are 



Energies and angular yields at the different ports from different PXR planes are shown in 
Tabled 

Table [6] shows the spectral brilliance expected from ports 4 and 5, the ones correspond¬ 
ing to the smallest Bragg angles and highest yields. With the crystal thickness kept constant 
at 0.168 mm, the spectral brilliance is highest for the (111) plane at these ports. The X-ray 
energy is between 7-8 keV and the ratio of the crystal thickness to attenuation length is 
close to the optimal value of around 0.2, seen in Figure [51 For the higher order planes, the 
PXR energy increases but the spectral brill lance decreases. Especially for planes (3,1,1) 
and (4,0,0) the brilliance drops by an order of magnitude compared to the (1,1,1) plane. 
This is partly due to the small relative thickness and increasing the crystal thickness would 
also increase the brilliance, but not significantly. The choice of plane would then be de¬ 
termined by whether higher energy or higher brilliance is more desirable. A higher energy 
beam, e.g. lOOMeV, would increase the yield and also the brilliance because the emittance 
growth due to multiple scattering would also be smaller. 
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Port 

Plane 111 

Plane 220 


Energy 

Ang. Yield 

Energy 

Ang. Yield 


[keV] 

[x 10^^ phot/e^-sr ] 

[keV] 

[x 10-5 phot/e--sr ] 

1 

4.3 

0.037 

6.9 

9.9 

2 

3.3 

3.4x10-5 

5.4 

1.2 

3 

3.2 

1.5x10-5 

5.3 

0.096 

4 

7.1 

76.9 

11.6 

90.5 

5 

7.9 

142 

12.8 

108.3 

Port 

P 

ane 311 

P 

ane 400 


Energy [keV] 

Ang. Yield 

Energy [keV] 

Ang. Yield 


[keV] 

[x 10-5 phot/e--sr ] 

[keV] 

[x 10-5 phot/e--sr ] 

1 

8.1 

5.3 

9.8 

8.8 

2 

6.4 

1.6 

7.7 

4.8 

3 

6.2 

1.4 

7.5 

4.6 

4 

13.6 

26.5 

16.4 

34.7 

5 

15.1 

41.3 

18.2 

40.8 


Table 5: X-ray energies at the different ports for different PXR reflection planes. Bragg 
geometry is assumed in all cases. 


Port 

Plane 111 

Plane 220 


Energy 

t ! La 


Energy 

t ! La 

Sp. Br. [xlO""] 


7.1 


5.5 

11.6 

0.09 

2.2 


7.9 

m 

5.7 

12.8 

0.07 

2.1 

Port 

Plane 311 

Plane 400 


Energy 

t / La 


Energy 

t / La 

Sp. Br. [xlO""] 



■oig 

0.59 

16.4 

0.04 

0.59 


■n 


0.57 

18.2 

0.03 

0.56 


Table 6: PXR photon energy, crystal thickness t relative to attenuation length La, and spec¬ 
tral brilliance (Sp. Br.) in units of photons/s-(mm-mrad)^- 0.1% BW] at ports 4 and 5 for 
different planes. The crystal thickness was 0.168mm in each case. 
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(110) plane 



Figure 7: Sketch of the (110) channeling plane, the (1,1,1) plane which is one of the 
possible PXR planes, the electron beam direction and the direction of the PXR photons 
from this reflection plane. 

6.2 PXR while channeling 

It was pointed out [|2Sl that if a beam is channeled within a crystal and emits channeling 
radiation, it may also emit PXR emission from reflection off complementary planes which 
intersect the channeling planes. It has subsequently been observed at the SAGA light source 
linac with 255 MeV electron beams [|29]|. Here we consider the prospect of detecting PXR 
emission under channeling conditions while using the present goniometer. As mentioned 
previously, this goniometer has a second port at 90° to the beam axis which could be used 
for detection of PXR. These requirements impose constraints on the possible PXR planes 
and the orientation of the crystal which we now consider. 

Choosing the (110) plane as the channeling plane, the unit normal to this plane is 
m = (1,1,0)/\/2. This vector must be normal to the velocity vector which could be of 
the form v = (a, —a,c)/\/2a^ + where (a,c) are arbitrary real numbers. We will con¬ 
sider two choices below, one with (a,c) = (—1,0) and another with (a,c) = (0,1). If the 
PXR plane is (111), the unit normal to the PXR plane isg = (1, —1,1) / \/3 and the Bragg 
angle 9b is determined by the condition sinOg = |g ■ v| = | (2a -|-c)/( ^^3{2a^ + c'^)) \. The 
choice V = (— l,l,0)/\/2 yields 9b = dsc^m{^j2/3) = 54.7° while choosing the vector 
V = (0,0,1) yields 9b = arcsin(A/l/3) = 35.3°. A smaller Bragg angle leads to a higher 
PXR photon energy, so we choose v = (0,0,1). In practice, with a given electron beam 
direction, these choices will correspond to different orientations of the crystal with respect 
to the beam velocity. For particles at other angles in the beam distribution, the requirement 
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for channeling is that the veloeity veetor has an angle smaller than the eritieal angle Oc, i.e. 
6 = aresin(m ■ v) < Oc- 

If we assume that the crystal has been eut so that the surfaee is parallel to the channeling 
plane, then the unit normal to the surface, defined as h (see Fig. O is identical to the veetor 
m defined above, i.e. n = m= (1,1,0)/ y/2. 


6.2.1 Crystal rotation 


If the erystal is aligned with so that the ehanneling plane is parallel to the veloeity veetor, 
then the Bragg angle with a ehosen PXR plane may not be appropriate to observe the 
PXR photons at 90 degrees to the beam. However, a rotation of the erystal about an axis 
orthogonal to the channeling plane can create the desired Bragg angle while maintaining 
the ehanneling eondition. If the rotation matrix about the normal to the ehanneling plane m 
is written as R(v/) where y/ is the rotation angle, then the rotated normal to the PXR plane 
is gR = R{yr)g. The requirement on this angle t// is that the velocity has the desired Bragg 
angle with the new normal gR, i.e. 


sin 05 = !/;•(-v) =R(v/)|-(-v) 


(28) 


The rotation matrix about an arbitrarily ehosen unit veetor fl = (jUx, Mz) is the so-ealled 
Rodrigues matrix lf30l given by 


R{xff) 


cosyr + n^S - sin yr + fly sin yr + 

Hz sin y/+ Hxl^yS cosy/+ ii^S -Hxsiny/+ HyHzS 

-jiy sin y/ + HxHzS Hx sin Y + ^y^z^ cos t// + ju/5 


(29) 


where S = 2sin^(t///2). When the ehanneling plane is (110), the rotation axis is fi = 
(1,1,0)/\/2. The desired angle of rotation Y is found from Eq. (l28l) with Qr = njA. 

Applying this to different possible PXR planes will ensure that the rotated light veetor 
lies in the (x,y) plane. A further rotation about the z axis (i.e along the beam direction) may 
be necessary to direct the light out along the x axis. When the emitted light as viewed at an 
angle 20g to the beam, the refleeted light veetor prior to any rotations is ^ = v — 2(g.v)g 
while that after rotation is Qr = v — 2(gR ■v)gR. 

As an example, eonsider the plane (1,1,1) for PXR production. Choosing the normal 
to this plane as g = (1, —1, —1)/\/3, we have 


gR ■ 


?) = ^(eos v/+\/2sinv/) = yr = 9JA°,CIr = ^(1,-1,0) 

vS \/2 \/2 


In this ease, a rotation by Y = 9.74° about the normal to the channeling plane and a further 
rotation by 0 = 45° about the beam axis will direct the PXR light out of the port along the 
positives axis. 

Table |7] shows some of the low order planes, the angles of rotation y and the 

energies of the emitted X-rays into the deteetor. This table shows that the (2,0,2) plane 
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Plane 

Angle v/[deg] 

Angle ^[deg] 

ExlkeV] 

(1,1,1) 

9.74 

45 

4.26 

(i,i,i) 

-9.74 

135 

4.26 

(2, 0. 2) 

0 

180 

6.95 

(2,0,2) 

0 

0 

6.95 

(0,2,2) 

0 

-90 

6.95 

(0,0,4) 

45 

135 

9.83 


Table 7: PXR planes, the angles of rotation y/ about the normal to the channeling plane 
(1,1,0) and (j), angle of rotation about the beam direction so that the PXR light is directed 
out of the port along the positive x direction at 90 degrees to the beam axis. The energy 
of the PXR light is also shown. Other low-order planes not shown are those for which no 
rotation will direct the X-ray beam in the desired direction. 

would be the simplest as it does not require any rotation of the crystal while it is oriented 
for channeling. Another reason for choosing this plane is that with rotation, the path length 
of the beam while channeling will be different compared to the unrotated case and could 
affect the channeling yield. 

The angle ^ between the unit normal h to the surface and the unit normal g to the 
channeling plane will change depending on the channeling plane chosen. We assume here 
as above that the crystal is cut parallel to the channeling plane so that n = (1,1,0) /\/2. The 
angle between the two is given by ^ = arccos(n ■ g), thus for the (2,0,2) PXR plane, this 
angle is 60°. 

6.2.2 Yield and spectral brilliance 

The above analysis has shown that with the crystal oriented for channeling along the (1,1,0) 
plane, PXR emission from the (2,0,2) plane can be obtained at 90° to the beam direction 
without any crystal rotation. In the subsequent discussion, we will primarily consider PXR 
from this plane. 

When beam energies are under 100 MeV, the channeling radiation spectrum consists of 
a few well defined lines and is best understood as arising from transitions between bound 
states of the transverse potential. We consider planar channeling and the x direction to be 
orthogonal to the channeling plane. The quantum mechanical states are then found from 
solutions of the one dimensional Schrodinger equation, 

+ = (30) 

Here V (x) is the one dimensional continuum potential obtained by averaging the three 
dimensional atomic potential along the orthogonal directions (y,z). 

Parametric X-rays emitted under channeling conditions (PXRC) are considered to be 
those in which the electrons stay within the same transverse energy band, i.e with only 
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intra-band transitions. PXR emitted while eleetrons transition between different energy 
bands is labeled as diffraeted ehanneling radiation (DCR) and has apparently not yet been 
experimentally observed. We will only eonsider the yield from PXRC here. 

The angular speetrum of PXRC is related to that of PXR by |[3T]| 


(fN ^ 


(31) 


where the sum over the states n ranges over the total number of bound states Ng- Here we 
also include the effects of a finite beam divergence, so we define the initial probability 
of occupation of state n by averaging over the beam divergence. It is given by 


{Pn) = 


1 1 
Vlnaidp 


J0;,,oexp[- 




x,0 


rdpl2 


2 ( 0 -;) 2 J-dp/2 


exp[-ik^j,^ox\ \l/n,Kc{^)dx\^)Kc (32) 


Here is the beam divergence in the channeling plane, dp is the inter-planar separation, k 
is the initial momentum wavenumber of the incident particle, ^x,o is the angle of incidence 
with respect to the channeling plane and ^fn,Kc is the wave function in the nth state with 
transverse wavenumber where —g/2 <Kc<g 2 - Since this wave function also depends 
upon the band wavenumber Kc in the Brillouin zone, the average on the right hand 
side represents an average over the wavenumbers in the Brillouin zone. 

The form factor describes the impact of the channeling wave functions on the PXR 
yield and is defined as Il28ll 

Fnn — { \ CXp [ ikxX\ \ ^fn) 

where k^ is the transverse component of the photon wave vector. We can therefore define 
the averaged form factor squared as 


{\Fnn{Ox)\^) = ^{\ [ V/„V(x)exp[-/— 


dl J-dp/2 


]yirn,Kcix)dx\ )kc 


(33) 


Here Ox is the angle of photon emission in the horizontal plane with respect to the direction 
of specular reflection and we average over the Brillouin zone as before. 


The quantum mechanical calculations were done with a Mathematica notebook used at 
the ELBE facility to study channeling and significantly corrected and modified for use 
at EAST, as described in [[3]|. The transverse function is expanded in a Eourier series using 
the lattice periodicity as 

oo 

Vi^)= Y, ^«exp[mgx] (34) 


The Eourier coefficients Vn are obtained from the Doyle-Tumer coefficients ll33ll 

as 


V„ = 


In 7,e 




«0 


\DT 


i=\ 


16%'^ 


{ngf 


(35) 


Here Vc is the volume of the unit cell, uq is the Bohr radius, Vj are the coordinates of the 
yth atom in the unit cell and M[g) = ^g^{u^) is the Debye-Waller factor (mentioned earlier 
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Figure 8: The differential yield as a funetion of the photon emission angle Ox (left) and 
Oy (right) without ehanneling and with ehanneling for three values of the initial beam di- 
vergenee. The ehanneling plane was (1,1,0) and the PXR plane was (2,0,2). Convolution 
broadening was ineluded in all eases. 


with Eq.© in Section |2l) describing thermal vibrations with mean squared amplitude (u^), 
assumed to be the same for all atoms. The wave functions are found by first expanding them 
in a series of Bloch functions and then solving the resulting matrix eigenvalue problem, see 

e.g [El- 

The wave functions are then subsequently used in calculating the correction factors 
defined above. The form factor (|E„„(0)p) = 1 and we find that in the range —S/y < 
Ox < 3 / 7 , (|E„„(0x)P) ~ 1 to within 0.2% in all cases. The significant correction for the 
PXRC yield is due to the initial probability of occupation (P„) which is determined by the 
beam divergence. This population of the bound states decreases as the beam divergence 
approaches the critical angle of channeling which at the FAST energy of 50MeV is about 
0.98mrad. The correction 5 is defined as the relative difference between the PXR and 
PXRC yields, hence as 5 = 1 — {Pn){\Fnn\'^)■ This correction factor increases with 
beam divergence and has the values 0.02, 0.09 and 0.58 at beam divergences of O.Olmrad, 
O.lmrad and Imrad respectively. 

Figure [8] shows the angular spectra from PXRC for the three divergences mentioned and 
compared with the spectrum from PXR without channeling. Table[8]shows that the angular 
yield decreases with beam divergence, because the fraction of particles in the bound states 
decreases with increasing divergence. Since the spectral brilliance is inversely proportional 
to the square of the spot size and assuming the emittance is conserved, a larger beam 
divergence implies a smaller spot size. The increase in brilliance with the smaller spot size 
dominates the decrease due to channeling. 

6.3 Possible applications of PXR 

One of the primary foreseen uses of PXR is in phase contrast imaging (PCI) of low Z 
materials, specially biological samples. There are many PCI methods, we plan to use the 
free space propagation method which does not need any additional optical elements. The 
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Plane 

Divergence 

AE 


Yield 

Spectral Brilliance 


mrad 

eV 

xl0~ 

^ photons/(e^-sr) 

photons/(mm-mrad)^-0.1% BW 




PXR 

PXRC 

PXR 

PXRC 

(2,0,2) 

0.01 

93 

0.40 

0.38 

35.9 

35.1 


0.1 

93 

0.40 

0.35 

2 .1x10^ 

1.9x10^ 


1.0 

93 

0.40 

0.16 

6.5x10^ 

1 .8x10^ 

(0,0,4) 

0.01 

131 

0.89 

0.85 

21.3 

20.8 


0.1 

131 

0.89 

0.78 

1.27 xlO^ 

1.14 xlO^ 


1.0 

131 


0.37 

4.25 xlO^ 

1 .2x10^ 


Table 8: Comparison of angular yields, photon flux and speetral brillianee without and with 
ehanneling for three values of the initial beam divergenee and for two PXR planes. The 
ehanneling plane was fixed at (1,1,0). Angle ^ = 60° for the (2,0,2) plane, and ^ = 90° 
(Laue geometry) for the (0,0,4) plane. The energy width is affeeted by the PXR plane but 
not by the initial beam divergenee. 


prineiple of PCI is that phase shifts undergone by hard X-rays when traversing low Z sam¬ 
ples are orders of magnitude larger than absorption effeets. The phase shift changes the 
complex amplitude of the wave field and hence causes intensity modulations when under¬ 
going interference with X-rays that have not passed through the sample. The quality of the 
image with the free space propagation method depends on the source size, the geometric 
magnification and the resolution of the detector. Th special feature of FAST that makes it 
suitable for PCI is the very low emittance of the electron beam, and hence the small X-ray 
source size after the crystal. With electron emittances around lOOnm using a conventional 
photocathode, beam sizes at the crystal around 1-5 jim can be achieved [[35l . This com¬ 
pares very favorably with recent PCI experiments using PXR lf36l which had beam sizes 
of 0.5 - 2 mm (FWHM) at the target. With the use of field emitter nanotips as the cathode, 
the beam emittance could be improved another order of magnitude IfTH implying a further 
improvement in image resolution. The use of the new goniometer, discussed in Section [Q] 
will enable imaging at multiple X-ray energies. 

Finally, we mention a proposal [[37l to generate short electron bunches at FAST using a 
slit mask placed in the middle of the bunch compressor chicane. Sub-picosecond electron 
bunches could be produced without the need of an undulator or an additional complex laser 
system. Due to scattering in the mask, the final beam intensity will only be about 10% 
of the initial intensity. However, starting from initial intensities of 1-3 nC, the final bunch 
intensities within the bunch train will be low enough (about 20 pC) for PXR generation. 
The resulting sub-picosecond X-ray pulses will have a higher peak brilliance and could be 
used for time resolved X-ray studies in materials science, chemistry and biology. 
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7 Conclusions 


In this paper we considered the prospect of generating PXR using diamond crystals with 50 
MeV electron beams from the photoinjector at the FAST facility at Fermilab. We revisited 
calculations of the energy width from both geometric and multiple scattering. Compar¬ 
isons with earlier experiments were found to yield reasonable agreement The PXR spec¬ 
trum model calculation was applied to the conditions at FAST. Using the presently available 
goniometer restricts the Bragg angle to 45° but allows a clear separation from the electron 
beam. The PXR energies from the planes studied fall in the range 4-10 keV with a spectral 
energy widths of ~ 1%. With a diamond crystal thickness of 168/im, maximum angular 
yields are about 10 ^ photons/e^-sr, taking into account attenuation in air from the crystal 
to the detector. Using crystals of different thicknesses, the PXR yield was found to saturate 
at a thickness of about l.2La, the photon attenuation length at two different energies. The 
spectral brilliance on the other hand attained a maximum value at around O.lLa for these 
same energies. This is mainly because the emittance growth over larger thicknesses reduces 
the brilliance more than the yield increases it. Next, PXR emission with a new goniometer 
with five possible ports was studied. The range of energies now spans 3-18 keV. The spec¬ 
tral brilliance, around 10^ photons/(mm-mrad)^-0.1% BW, is reached with a fixed crystal 
thickness which is less than O.lLa for most energies, so higher brilliance is feasible with 
thicker crystals. Use of this goniometer would open up the possibility of extracting PXR 
at multiple energies simultaneously. Finally, PXRC or PXR under channeling conditions 
was studied and the yield with quantum corrections from channeling was calculated for 
three beam divergences. While the reduction in PXR yield during channeling was smallest 
for the lowest divergence, higher brilliance favors the smallest beam spot size or equiva¬ 
lently the largest divergence under conditions of equal emittances. This PXRC emission 
makes possible simultaneous X-ray emission from channeling and PXR at 90 degrees to 
each other. The brilliance of PXR appears to be sufficient for phase contrast imaging and 
the FAST facility with (10 - 100) nm scale electron emittances should enable imaging with 
very good resolution. 


Acknowledgments 

We thank the Lee Teng undergraduate internship program at Fermilab which awarded T. 
Seiss a summer internship in 2014. Fermilab is operated by Fermi Research Alliance under 
DOE Contract No. DE-AC02-07-CH11359. 


References 

[1] E. Harms et al, ICEA Beam Dynamics Newslett. 64, 133 (2014); 

[2] D. Mihalcea et al, Proc. IPAC15 (2015) 

[3] T. Sen and C. Eynn, Int. J. Mod. Phys. A, 29, 1450179 (2014) 


23 



[4] M.L. Ter-Mikaelian, High energy electromagnetic processes in condensed media, Wi¬ 
ley Interscience, New York (1972) 

[5] V.G. Baryshevsky and I.D. Feranchuk, Sov. Phys. JETP, 34, 502 (1972) 

[6] A.N. Didenko et al, Phys. Lett. A, 110, 177 (1985) 

[7] P. Rullhusen et al. Novel Radiation Sources using Relativistic Electrons, World Scien¬ 
tific Publishing Co. (1998) 

[8] V.G. Baryshevsky et al. Parametric X-Ray Radiation in Crystals, Springer (2005) 

[9] W. Scandale et al, Phys. Lett. B, 701, 180 (2011) 

[10] W. Gabella et al, Nucl. Instrum, & Meth. B309, 10 (2013) 

[11] P. Plot et al, Appl. Phys. Lett., 104, 263504 (2014) 

[12] I.D. Leranchuk and I.V. Ivashin, J. Phys (Paris), 46, 1981 (1985) 

[13] H. Nitta, Phys. Lett. A, 158, 270 (1991) 

[14] K. Brenzinger et al, Z. Phys. A,358, 107 (1997) 

[15] D.T. Cromer and J.B. Mann. Act Crys. A24, 321 (1968) 

[16] C.T. Chantler, X-ray form factor, attenuation and scatter¬ 
ing tables,J. Phys.Chem.Ref.Data 24,71 (1995); Table at http: 

physics.nist.gov/PhysRefData/PPast/html/form.html 

[17] V.G. Baryshevsky and I.D. Leranchuk, Sov. Phys. JETP 34, 502 (1972) 

[18] G.M. Garibian and C. Yang, Sov. Phys. JETP, 34, 495 (1972) 

[19] A. Caticha, Phys. Rev. A, 40, 4322 (1989) 

[20] H. Backe et al, Nucl. Instr. & Meth. B, 234, 138 (2005) 

[21] J. Beringer et al (Particle Data Group), Phys. Rev. D,86, 010001 (2012) 

[22] A. Caticha, Phys. Rev. B, 45, 9541 (1992) 

[23] J. Preudenberger et al, Appl. Phys. Lett. 70,267 (1997) 

[24] B Sones PhD thesis, Rensselaer Polytechnic Inst., Troy, NY (2004) 

[25] B. Sones et al, Nucl. Instr & Meth A, 560, 589 (2006) 

[26] W. Wagner et al, Nucl. Instr. & Meth. B266, 327 (2008) 

[27] R.K. Li et al, Phys. Rev. ST-AB, 15, 090702 (2012) 

[28] R. Yabuki et al, Phys. Rev. B, 63, 174112 (2001) 


24 



[29] Y. Takabayashi et al, Nucl. Instr. & Meth. B, 315, 105 (2013) 

[30] O. Rodrigues, J. Math. Pures. Appl. 5, 380 (1840) 

[31] K.B. Korotchenko et al, JETP, 95, 481 (2012) 

[32] B. Azadegan, Comp. Phys. Comm., 184, 1064 (2013) 

[33] R.A. Doyle and PS. Turner, Aeta Crystallogr. A , 24, 390 (1968) 

[34] B. Azadegan et al, Phys. Rev B.,74, 045209 (2006) 

[35] P. Piot et al, AIP Conf. Proe. 1507, 732 (2012) 

[36] Y. Hayakawa et al, J. Inst. 8, C08001 (2013) 

[37] J.C. Thangaraj and P. Piot, arXiv:1310.5389 


25 


